********************************************************************************
* 								 MAIN ANALYSES								   *
********************************************************************************

clear all
set more off
set scheme s2color
grstyle init
grstyle set plain, horizontal grid dotted box
macro drop _all
est drop _all
set matsize 800
set seed 8675309

** Set Working Directory

if c(username) == "cb2257" {
global dir "~/Desktop"  // CWB's Directory
global main "${dir}/APSR Replication"
global code "${main}/Analysis"	
global data "${main}/Data"	
global figs "${main}/Figures"
global tabs "${main}/Tables"																			
}

if c(username) == "christopherblair" {
global dir "~/Desktop"  // CWB's Directory
global main "${dir}/APSR Replication"
global code "${main}/Analysis"	
global data "${main}/Data"	
global figs "${main}/Figures"
global tabs "${main}/Tables"																			
}

else if c(username) == "yourusername"{
global dir "~/yourfilepath"  // YOUR Directory
}

cd "${main}"

********************************************************************************
********************************************************************************
********************************************************************************
********************************************************************************
********************************************************************************

use "${data}/Survey/survey_final.dta", clear
est drop _all

sort respondent_id

global core "republican democrat woman highschool college age"
global move "owner timecomm community_links"

qui reghdfe personal_hurricane 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
gen regsample=(e(sample))

** Table SI-1 of the Additional Supplement
estpost sum race_white race_black race_latinx race_multi race_asian race_indigenous woman age highschool college income republican democrat if regsample==1
esttab using "${tabs}/SumStat.tex", style(tex) cells("count(pattern(1 1 0) fmt(0)) mean(pattern(1 1 0) fmt(3)) sd(pattern(1 1 0) fmt(3)) min(pattern(1 1 0) fmt(3)) max(pattern(1 1 0) fmt(3))") varlabels(race_white "& \\ \textsc{Demographics} \\ \hspace{3mm}Race: White" race_black "\hspace{3mm}Race: Black" race_latinx "\hspace{3mm}Race: Latinx" race_multi "\hspace{3mm}Race: Multiracial" race_asian "\hspace{3mm}Race: Asian" race_indigenous "\hspace{3mm}Race: Indigenous" woman "\hspace{3mm}Woman" age "\hspace{3mm}Age" highschool "\hspace{3mm}Education: High Scool Graduate" college "\hspace{3mm}Education: College Graduate" income "\hspace{3mm}Income" republican "\hspace{3mm}Party: Republican" democrat "\hspace{3mm}Party: Democrat") nolines prehead(\begin{tabular}{l*{7}{c}} \hline & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) collabels("Observations" "Mean" "Std. Dev." "Min" "Max") nomti nonum noobs replace

estpost sum race_white race_black race_latinx race_multi race_asian race_indigenous woman age highschool college income republican democrat if regsample==1 [w=eweight]
esttab using "${tabs}/SumStat_Nation.tex", style(tex) cells("count(pattern(1 1 0) fmt(0)) mean(pattern(1 1 0) fmt(3)) sd(pattern(1 1 0) fmt(3)) min(pattern(1 1 0) fmt(3)) max(pattern(1 1 0) fmt(3))") varlabels(race_white "& \\ \textsc{Demographics} \\ \hspace{3mm}Race: White" race_black "\hspace{3mm}Race: Black" race_latinx "\hspace{3mm}Race: Latinx" race_multi "\hspace{3mm}Race: Multiracial" race_asian "\hspace{3mm}Race: Asian" race_indigenous "\hspace{3mm}Race: Indigenous" woman "\hspace{3mm}Woman" age "\hspace{3mm}Age" highschool "\hspace{3mm}Education: High Scool Graduate" college "\hspace{3mm}Education: College Graduate" income "\hspace{3mm}Income" republican "\hspace{3mm}Party: Republican" democrat "\hspace{3mm}Party: Democrat") nolines prehead(\begin{tabular}{l*{7}{c}} \hline & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) collabels("Observations" "Mean" "Std. Dev." "Min" "Max") nomti nonum noobs replace

estpost sum race_white race_black race_latinx race_multi race_asian race_indigenous woman age highschool college income republican democrat if regsample==1 [w=rweight]
esttab using "${tabs}/SumStat_Region.tex", style(tex) cells("count(pattern(1 1 0) fmt(0)) mean(pattern(1 1 0) fmt(3)) sd(pattern(1 1 0) fmt(3)) min(pattern(1 1 0) fmt(3)) max(pattern(1 1 0) fmt(3))") varlabels(race_white "& \\ \textsc{Demographics} \\ \hspace{3mm}Race: White" race_black "\hspace{3mm}Race: Black" race_latinx "\hspace{3mm}Race: Latinx" race_multi "\hspace{3mm}Race: Multiracial" race_asian "\hspace{3mm}Race: Asian" race_indigenous "\hspace{3mm}Race: Indigenous" woman "\hspace{3mm}Woman" age "\hspace{3mm}Age" highschool "\hspace{3mm}Education: High Scool Graduate" college "\hspace{3mm}Education: College Graduate" income "\hspace{3mm}Income" republican "\hspace{3mm}Party: Republican" democrat "\hspace{3mm}Party: Democrat") nolines prehead(\begin{tabular}{l*{7}{c}} \hline & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) collabels("Observations" "Mean" "Std. Dev." "Min" "Max") nomti nonum noobs replace

drop regsample
eststo clear


********************************************************************************
* 							 EXPERIENCED HURRICANE							   *
********************************************************************************
** Table SI-2 of the Additional Supplement

eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe personal_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=ian_windswath2
eststo: reghdfe personal_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Windswath"

replace exposure=i_ian_surge
eststo: reghdfe personal_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Storm Surge"

replace exposure=fema_aid
eststo: reghdfe personal_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "FEMA Aid"

replace exposure=hurricane_index1
eststo: reghdfe community_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=ian_windswath2
eststo: reghdfe community_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Windswath"

replace exposure=i_ian_surge
eststo: reghdfe community_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Storm Surge"

replace exposure=fema_aid
eststo: reghdfe community_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "FEMA Aid"

esttab using "${tabs}/Experience_Hurricane.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{Personally Experienced a}\\\textbf{Hurricane in Past Year (=1)}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Community Experienced a}\\\textbf{Hurricane in Past Year (=1)}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 5) replace
esttab using "${tabs}/Experience_Hurricane_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{Personally Experienced a}\\\textbf{Hurricane in Past Year (=1)}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Community Experienced a}\\\textbf{Hurricane in Past Year (=1)}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 5) replace

eststo clear

*** OTHER FORMS OF CLIMATE CHANGE? ***
** Table SI-3 of the Additional Supplement

eststo clear

replace exposure=hurricane_index1
eststo: reghdfe personal_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe personal_flood 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe personal_fire 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe personal_drought 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

esttab using "${tabs}/Experience_OtherClimate.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\shortstack{\textbf{Personally Experienced Extreme Weather in Past Year (=1)}}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Hurricane" "Floods" "Wildfires" "Drought") extracols(1 2 3 4) replace
esttab using "${tabs}/Experience_OtherClimate_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\shortstack{\textbf{Personally Experienced Extreme}\\\textbf{Weather in Past Year (=1)}}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Hurricane" "Floods" "Wildfires" "Drought") extracols(1 2 3 4) replace

drop exposure
eststo clear

** Figure 3 of the Manuscript
eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe personal_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=hurricane_index1
eststo: reghdfe community_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe personal_flood 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe personal_fire 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe personal_drought 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est2, msymbol(D) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est1, msymbol(O) mfcolor(black%0) mlcolor(black%0) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black%0) mlabpos(1) mlabgap(*2) ciopts(lcolor(black%0 black%0) lwidth(.55 1.1))) (est3, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(2) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est4, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est5, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.65 `" "Personally Experienced a" "Hurricane in Past Year (=1)" "' .79 `" "Community Experienced a" "Hurricane in Past Year (=1)" "' 1.07 `" "Personally Experienced a" "Flood in Past Year (=1)" "' 1.21 `" "Personally Experienced a" "Wildfire in Past Year (=1)" "' 1.36 `" "Personally Experienced a" "Drought in Past Year (=1)" "') xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xlabel(-.05(.05).25) xmtick(-.05(.025).25) legend(off) yline(.9, lcolor(black) lwidth(*1.5)) xline(0, lcolor(cranberry) lpatt(shortdash))
graph export "${figs}/experience_weather.png", replace

drop exposure
eststo clear

********************************************************************************
* 					 EXPERIENCED CLIMATE MIGRATION							   *
********************************************************************************
** Table SI-4 of the Additional Supplement

eststo clear

gen hurricane_personal=climmig_personal*climmig_reason_hurricane
gen hurricane_family=climmig_family*climmig_reason_hurricane
gen hurricane_friend=climmig_friend*climmig_reason_hurricane
gen hurricane_acquaintance=climmig_acquaintance*climmig_reason_hurricane

gen exposure=hurricane_index1 
eststo: reghdfe climmig_reason_hurricane 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=hurricane_index1
eststo: reghdfe climmig_reason_flood 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=hurricane_index1
eststo: reghdfe climmig_reason_fire 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=hurricane_index1
eststo: reghdfe climmig_reason_drought 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=hurricane_index1
eststo: reghdfe hurricane_personal 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=hurricane_index1
eststo: reghdfe hurricane_family 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=hurricane_index1
eststo: reghdfe hurricane_friend 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=hurricane_index1
eststo: reghdfe hurricane_acquaintance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

esttab using "${tabs}/Experience_ClimateMigration.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{I or Someone I Know Moved}\\\textbf{for Climate-Related Reasons (=1)}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Know Someone Who Moved}\\\textbf{Because of a Hurricane (=1)}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Hurricane" "Floods" "Wildfires" "Drought" "Personally Moved" "Family Moved" "Friend Moved" "Acquaintance Moved") extracols(1 5) replace
esttab using "${tabs}/Experience_ClimateMigration_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{I or Someone I Know Moved}\\\textbf{for Climate-Related Reasons (=1)}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Know Someone Who Moved}\\\textbf{Because of a Hurricane (=1)}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Hurricane" "Floods" "Wildfires" "Drought" "Personally Moved" "Family Moved" "Friend Moved" "Acquaintance Moved") extracols(1 5) replace

drop exposure
eststo clear

********************************************************************************
* 						 PLACEBO: MILITARY SUPPORT							   *
********************************************************************************
** Figure A-10 of the Appendix and Table D-30 of the Additional Supplement

eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe priority_military_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Exposure "Index"

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(1 `" "Strengthening the Military" "is a Top Priority (=1)" "') xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xlabel(-.04(.02).08) xmtick(-.04(.01).08) legend(off) xline(0, lcolor(cranberry) lpatt(shortdash)) ysize(1) xsize(1.75)
graph export "${figs}/placebo_outcome.png", replace

eststo: reghdfe priority_economy_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Exposure "Index"

eststo: reghdfe priority_migration_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Exposure "Index"

eststo: reghdfe priority_health_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Exposure "Index"

esttab using "${tabs}/Placebo_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{6}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{Perceived Issue Importance of:}}} \\ \cline{3-6} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{Strengthening\\Military}" "\shortstack{Strengthening\\Economy}" "\shortstack{Addressing\\Migration}" "\shortstack{Strengthening\\Healthcare}") extracols(1) replace
drop exposure
eststo clear

********************************************************************************
*		 							 MAIN OUTCOMES							   *
********************************************************************************
** Table 2 of the Main Text

eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe climmig_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climmig_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_mitigation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_adaptation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_science 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

esttab using "${tabs}/Climate_Outcomes_Main.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Migration}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}\\\textbf{Policies}}} & & \shortstack{\textbf{Science of}\\\textbf{Climate Change}}} \\ \cline{3-4} \cline{6-7} \cline{9-10} \cline{12-12} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "Mitigation" "Adaptation" "Science") extracols(1 3 5 7) replace
esttab using "${tabs}/Climate_Outcomes_Main_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "& \\ Democrat" woman "& \\ Woman" highschool "& \\ High School Graduate" college "& \\ College Graduate" age "& \\ Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Migration}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}\\\textbf{Policies}}} & & \shortstack{\textbf{Science of}\\\textbf{Climate Change}}} \\ \cline{3-4} \cline{6-7} \cline{9-10} \cline{12-12} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{Issue Importance}" "\shortstack{Policy Action}" "\shortstack{Issue Importance}" "\shortstack{Policy Action}" "Mitigation" "Adaptation" "Science") extracols(1 3 5 7) replace

matrix Means = J(7, 5, .) 
matrix list Means 
est resto est1
matrix Means[1, 1] = _b[1.post#c.exposure]
est resto est2
matrix Means[2, 1] = _b[1.post#c.exposure]
est resto est3
matrix Means[3, 1] = _b[1.post#c.exposure]
est resto est4
matrix Means[4, 1] = _b[1.post#c.exposure]
est resto est5
matrix Means[5, 1] = _b[1.post#c.exposure]
est resto est6
matrix Means[6, 1] = _b[1.post#c.exposure]
est resto est7
matrix Means[7, 1] = _b[1.post#c.exposure]
eststo clear
matrix Means[1, 3] = .0041827
matrix Means[2, 3] = -.0388997
matrix Means[3, 3] = .0110551
matrix Means[4, 3] = -.0077864
matrix Means[5, 3] = -.0548088
matrix Means[6, 3] = -.0529555
matrix Means[7, 3] = -.0104822

drop exposure
eststo clear

********************************************************************************
*		 							 OSTER'S BOUNDS							   *
********************************************************************************
** Table A-13 of the Main Appendix

eststo clear

gen exposure=hurricane_index1 
reghdfe climmig_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
gen estimation=e(sample)

matrix Oster = J(7, 1, .) 
matrix list Oster 

eststo: areg climmig_importance 1.post#c.exposure $core i.start_date if estimation==1 [pw=eweight], cluster(geoid) abs(geoid)
psacalc delta 1.post#c.exposure, rmax(.21632)
matrix Oster[1, 1] = `r(delta)'

eststo: areg climmig_action 1.post#c.exposure $core i.start_date if estimation==1 [pw=eweight], cluster(geoid) abs(geoid)
psacalc delta 1.post#c.exposure, rmax(.40157)
matrix Oster[2, 1] = `r(delta)'

eststo: areg climchg_importance 1.post#c.exposure $core i.start_date if estimation==1 [pw=eweight], cluster(geoid) abs(geoid)
psacalc delta 1.post#c.exposure, rmax(.3354)
matrix Oster[3, 1] = `r(delta)'

eststo: areg climchg_action 1.post#c.exposure $core i.start_date if estimation==1 [pw=eweight], cluster(geoid) abs(geoid)
psacalc delta 1.post#c.exposure, rmax(.36192)
matrix Oster[4, 1] = `r(delta)'

eststo: areg climchg_mitigation 1.post#c.exposure $core i.start_date if estimation==1 [pw=eweight], cluster(geoid) abs(geoid)
psacalc delta 1.post#c.exposure, rmax(.39598)
matrix Oster[5, 1] = `r(delta)'

eststo: areg climchg_adaptation 1.post#c.exposure $core i.start_date if estimation==1 [pw=eweight], cluster(geoid) abs(geoid)
psacalc delta 1.post#c.exposure, rmax(.30498)
matrix Oster[6, 1] = `r(delta)'

eststo: areg climchg_science 1.post#c.exposure $core i.start_date if estimation==1 [pw=eweight], cluster(geoid) abs(geoid)
psacalc delta 1.post#c.exposure, rmax(.29237)
matrix Oster[7, 1] = `r(delta)'

matrix list Oster

drop exposure 
eststo clear

********************************************************************************
*		 					 CINELLI AND HAZLETT BOUNDS						   *
********************************************************************************
** Figure A-7 of the Main Appendix

eststo clear

gen exposure=binary_index1*post
gen Democrat=democrat
global chcore "republican Democrat woman highschool college age"

sensemakr climmig_importance exposure $chcore i.geoid i.start_date if estimation==1, treat(exposure) benchmark(Democrat) contourplot clim(0 .161) kd(1 2 3) c(6)
graph save "${figs}/chbounds_climmig_importance.gph", replace
*graph export "${figs}/chbounds_climmig_importance.png", replace

sensemakr climmig_action exposure $chcore i.geoid i.start_date if estimation==1, treat(exposure) benchmark(Democrat) contourplot clim(0 .161) kd(1 2 3) c(6)
graph save "${figs}/chbounds_climmig_action.gph", replace
*graph export "${figs}/chbounds_climmig_action.png", replace

sensemakr climchg_importance exposure $chcore i.geoid i.start_date if estimation==1, treat(exposure) benchmark(Democrat) contourplot clim(0 .161) kd(1 2 3) c(6)
graph save "${figs}/chbounds_climchg_importance.gph", replace
*graph export "${figs}/chbounds_climchg_importance.png", replace

sensemakr climchg_action exposure $chcore i.geoid i.start_date if estimation==1, treat(exposure) benchmark(Democrat) contourplot clim(0 .161) kd(1 2 3) c(6)
graph save "${figs}/chbounds_climchg_action.gph", replace
*graph export "${figs}/chbounds_climchg_action.png", replace

sensemakr climchg_mitigation exposure $chcore i.geoid i.start_date if estimation==1, treat(exposure) benchmark(Democrat) contourplot clim(0 .161) kd(1 2 3) c(6)
graph save "${figs}/chbounds_climchg_mitigation.gph", replace
*graph export "${figs}/chbounds_climchg_mitigation.png", replace

sensemakr climchg_adaptation exposure $chcore i.geoid i.start_date if estimation==1, treat(exposure) benchmark(Democrat) contourplot clim(0 .161) kd(1 2 3) c(6)
graph save "${figs}/chbounds_climchg_adaptation.gph", replace
*graph export "${figs}/chbounds_climchg_adaptation.png", replace

sensemakr climchg_science exposure $chcore i.geoid i.start_date if estimation==1, treat(exposure) benchmark(Democrat) contourplot clim(0 .161) kd(1 2 3) c(6)
graph save "${figs}/chbounds_climchg_science.gph", replace
*graph export "${figs}/chbounds_climchg_science.png", replace

drop exposure estimation
eststo clear

********************************************************************************
*		 					 MAIN OUTCOMES -- NO FE							   *
********************************************************************************
** Result not presented in paper

eststo clear

reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
gen sample=(e(sample))

gen exposure=hurricane_index1 
eststo: reghdfe climmig_importance 1.post#c.exposure 1.post c.exposure $core if sample==1 [pw=eweight], cluster(geoid) noabs
estadd local County "No"
estadd local Date "No"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climmig_action 1.post#c.exposure 1.post c.exposure $core if sample==1 [pw=eweight], cluster(geoid) noabs
estadd local County "No"
estadd local Date "No"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_importance 1.post#c.exposure 1.post c.exposure $core if sample==1 [pw=eweight], cluster(geoid) noabs
estadd local County "No"
estadd local Date "No"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_action 1.post#c.exposure 1.post c.exposure $core if sample==1 [pw=eweight], cluster(geoid) noabs
estadd local County "No"
estadd local Date "No"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_mitigation 1.post#c.exposure 1.post c.exposure $core if sample==1 [pw=eweight], cluster(geoid) noabs
estadd local County "No"
estadd local Date "No"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_adaptation 1.post#c.exposure 1.post c.exposure $core if sample==1 [pw=eweight], cluster(geoid) noabs
estadd local County "No"
estadd local Date "No"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_science 1.post#c.exposure 1.post c.exposure $core if sample==1 [pw=eweight], cluster(geoid) noabs
estadd local County "No"
estadd local Date "No"
estadd local Demographic "Yes"
estadd local Exposure "Index"

esttab using "${tabs}/Climate_Outcomes_NoFE.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure 1.post exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" 1.post "& \\ Post" exposure "& \\ Hurricane Exposure") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Migration}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}\\\textbf{Policies}}} & & \shortstack{\textbf{Science of}\\\textbf{Climate Change}}} \\ \cline{3-4} \cline{6-7} \cline{9-10} \cline{12-12} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "Mitigation" "Adaptation" "Science") extracols(1 3 5 7) replace
esttab using "${tabs}/Climate_Outcomes_NoFE_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure 1.post exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" 1.post "& \\ Post" exposure "& \\ Hurricane Exposure" republican "& \\ Republican" democrat "& \\ Democrat" woman "& \\ Woman" highschool "& \\ High School Graduate" college "& \\ College Graduate" age "& \\ Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Migration}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}\\\textbf{Policies}}} & & \shortstack{\textbf{Science of}\\\textbf{Climate Change}}} \\ \cline{3-4} \cline{6-7} \cline{9-10} \cline{12-12} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{Issue Importance}" "\shortstack{Policy Action}" "\shortstack{Issue Importance}" "\shortstack{Policy Action}" "Mitigation" "Adaptation" "Science") extracols(1 3 5 7) replace

drop exposure
eststo clear

********************************************************************************
*		 							 PCF OUTCOMES							   *
********************************************************************************
** Table A-5 of the Main Appendix

eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe climmig_importance_pcf 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climmig_action_pcf 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_importance_pcf 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_action_pcf 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_mitigation_pcf 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_adaptation_pcf 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_science_pcf 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

esttab using "${tabs}/Climate_Outcomes_PCF.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Migration}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}\\\textbf{Policies}}} & & \shortstack{\textbf{Science of}\\\textbf{Climate Change}}} \\ \cline{3-4} \cline{6-7} \cline{9-10} \cline{12-12} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "Mitigation" "Adaptation" "Science") extracols(1 3 5 7) replace
esttab using "${tabs}/Climate_Outcomes_PCF_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Migration}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}\\\textbf{Policies}}} & & \shortstack{\textbf{Science of}\\\textbf{Climate Change}}} \\ \cline{3-4} \cline{6-7} \cline{9-10} \cline{12-12} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "Mitigation" "Adaptation" "Science") extracols(1 3 5 7) replace

drop exposure
eststo clear

********************************************************************************
*		 							 EW OUTCOMES							   *
********************************************************************************
** Table A-6 of the Main Appendix

eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe climmig_importance_ew 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climmig_action_ew 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_importance_ew 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_action_ew 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_mitigation_ew 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_adaptation_ew 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_science_ew 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

esttab using "${tabs}/Climate_Outcomes_EW.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Migration}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}\\\textbf{Policies}}} & & \shortstack{\textbf{Science of}\\\textbf{Climate Change}}} \\ \cline{3-4} \cline{6-7} \cline{9-10} \cline{12-12} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "Mitigation" "Adaptation" "Science") extracols(1 3 5 7) replace
esttab using "${tabs}/Climate_Outcomes_EW_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Migration}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}}} & & \multicolumn{2}{c}{\shortstack{\textbf{Climate Change}\\\textbf{Policies}}} & & \shortstack{\textbf{Science of}\\\textbf{Climate Change}}} \\ \cline{3-4} \cline{6-7} \cline{9-10} \cline{12-12} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "\shortstack{ \\Issue\\Importance}" "\shortstack{ \\Policy\\Action}" "Mitigation" "Adaptation" "Science") extracols(1 3 5 7) replace

drop exposure
eststo clear

********************************************************************************
* 							 CLIMATE MIGRATION								   *
********************************************************************************
** Table SI-5 of the Additional Supplement

eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe climmig_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=ian_windswath2
eststo: reghdfe climmig_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Windswath"

replace exposure=i_ian_surge
eststo: reghdfe climmig_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Storm Surge"

replace exposure=fema_aid
eststo: reghdfe climmig_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "FEMA Aid"

replace exposure=hurricane_index1
eststo: reghdfe climmig_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=ian_windswath2
eststo: reghdfe climmig_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Windswath"

replace exposure=i_ian_surge
eststo: reghdfe climmig_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Storm Surge"

replace exposure=fema_aid
eststo: reghdfe climmig_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "FEMA Aid"

esttab using "${tabs}/Climate_Migration.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{Issue Importance of Climate Migration}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Policy Action on Climate Migration}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 5) replace
esttab using "${tabs}/Climate_Migration_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{Issue Importance of Climate Migration}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Policy Action on Climate Migration}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 5) replace

drop exposure
eststo clear

** Figure 4a of the Main Paper
gen exposure=hurricane_index1 
eststo: reghdfe climmig_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe priority_climmig_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatemig1_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatemig2_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est1, msymbol(O) mfcolor(black%0) mlcolor(black%0) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black%0) mlabpos(1) mlabgap(*2) ciopts(lcolor(black%0 black%0) lwidth(.55 1.1))) (est2, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est3, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est4, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.67 `" "Issue Importance of Climate Migration" "' 1 `" "Addressing Climate Migration" "is a Top Priority (=1)" "' 1.17 `" "Climate Migration is" "a Serious Problem (=1)" "' 1.34 `" "Climate Migration Will Have a" "Serious Impact During My Life (=1)" "') xtitle("Hurricane Exposure x Post", size(vlarge)) xscale(titlegap(*5)) xlabel(-.04(.04).2) xmtick(-.04(.02).2) legend(off) yline(.835, lcolor(black) lwidth(*1.5)) xline(0, lcolor(cranberry) lpatt(shortdash))
graph export "${figs}/climmig_importance.png", replace
esttab using "${tabs}/ClimMig_Importance_SubItems.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\shortstack{\textbf{Issue Importance of Climate Migration}}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Index" "\shortstack{Addressing Climate Migration\\is a Top Priority}" "\shortstack{Climate Migration is\\a Serious Problem}" "\shortstack{Climate Migration Will\\Have a Serious Impact\\During my Life}") extracols(1 2 3 4) replace

eststo clear
** Figure 4b of the Main Paper

eststo: reghdfe climmig_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_usclimmig_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_intlclimmig_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatemig3_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatemig4_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatemig5_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est1, msymbol(O) mfcolor(black%0) mlcolor(black%0) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black%0) mlabpos(1) mlabgap(*2) ciopts(lcolor(black%0 black%0) lwidth(.55 1.1))) (est2, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est3, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est4, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est5, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est6, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.63 `" "Policy Action on Climate Migration" "' .88 `" "Tax Increase to Resettle" "Internal Climate Migrants (=1)" "' 1 `" "Tax Increase to Resettle" "International Climate Migrants (=1)" "' 1.13 `" "Vote for a Politician Who Promised" "to Address Climate Migration (=1)" "' 1.25 `" "U.S. Should Do More" "to Help Climate Migrants (=1)" "' 1.38 `" "International Community Should Do" "More to Help Climate Migrants (=1)" "') xtitle("Hurricane Exposure x Post", size(vlarge)) xscale(titlegap(*5)) xlabel(-.04(.04).2) xmtick(-.04(.02).2) legend(off) yline(.745, lcolor(black) lwidth(*1.5)) xline(0, lcolor(cranberry) lpatt(shortdash))
graph export "${figs}/climmig_action.png", replace
esttab using "${tabs}/ClimMig_Action_SubItems.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{13}{c}} \hline & \\ & & \multicolumn{11}{c}{\shortstack{\textbf{Policy Action on Climate Migration}}} \\ \cline{3-13} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Index" "\shortstack{Tax Increase to\\Resettle Internal\\Climate Migrants}" "\shortstack{Tax Increase to\\Resettle Int'l\\Climate Migrants}" "\shortstack{Vote for a Politician\\Who Promised to\\Address Climate Migration}" "\shortstack{U.S. Should Do\\More to Help\\Climate Migrants}" "\shortstack{Int'l Community Should\\Do More to Help\\Climate Migrants}") extracols(1 2 3 4 5 6) replace

drop exposure
eststo clear

********************************************************************************
* 								 CLIMATE CHANGE								   *
********************************************************************************
** Table SI-6 of the Additional Supplement

eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe climchg_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=ian_windswath2
eststo: reghdfe climchg_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Windswath"

replace exposure=i_ian_surge
eststo: reghdfe climchg_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Storm Surge"

replace exposure=fema_aid
eststo: reghdfe climchg_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "FEMA Aid"

replace exposure=hurricane_index1
eststo: reghdfe climchg_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=ian_windswath2
eststo: reghdfe climchg_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Windswath"

replace exposure=i_ian_surge
eststo: reghdfe climchg_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Storm Surge"

replace exposure=fema_aid
eststo: reghdfe climchg_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "FEMA Aid"

esttab using "${tabs}/Climate_Change.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{Issue Importance of Climate Change}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Policy Action on Climate Change}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 5) replace
esttab using "${tabs}/Climate_Change_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{Issue Importance of Climate Change}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Policy Action on Climate Change}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 5) replace

drop exposure
eststo clear
** Figure 4c of the Main Paper

gen exposure=hurricane_index1 
eststo: reghdfe climchg_importance 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe priority_climchg_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatechg1_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatechg2_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est1, msymbol(O) mfcolor(black%0) mlcolor(black%0) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black%0) mlabpos(1) mlabgap(*2) ciopts(lcolor(black%0 black%0) lwidth(.55 1.1))) (est2, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est3, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est4, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.67 `" "Issue Importance of Climate Change" "' 1 `" "Addressing Climate Change" "is a Top Priority (=1)" "' 1.17 `" "Climate Change is" "a Serious Problem (=1)" "' 1.34 `" "Climate Change Will Have a" "Serious Impact During My Life (=1)" "') xtitle("Hurricane Exposure x Post", size(vlarge)) xscale(titlegap(*5)) xlabel(-.04(.04).2) xmtick(-.04(.02).2) legend(off) yline(.835, lcolor(black) lwidth(*1.5)) xline(0, lcolor(cranberry) lpatt(shortdash))
graph export "${figs}/climchg_importance.png", replace
esttab using "${tabs}/ClimChg_Importance_SubItems.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\shortstack{\textbf{Issue Importance of Climate Change}}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Index" "\shortstack{Addressing Climate Change\\is a Top Priority}" "\shortstack{Climate Change is\\a Serious Problem}" "\shortstack{Climate Change Will\\Have a Serious Impact\\During my Life}") extracols(1 2 3 4) replace

eststo clear
** Figure 4d of the Main Paper

eststo: reghdfe climchg_action 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatechg3_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatechg4_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatechg5_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climatechg6_high 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est1, msymbol(O) mfcolor(black%0) mlcolor(black%0) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black%0) mlabpos(1) mlabgap(*2) ciopts(lcolor(black%0 black%0) lwidth(.55 1.1))) (est2, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est3, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est4, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est5, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.63 `" "Policy Action on Climate Change" "' .93 `" "Tax Increase to Fund Programs" "to Reduce Climate Change (=1)" "' 1.07 `" "Vote for a Politician Who Promised" "to Reduce Climate Change (=1)" "' 1.22 `" "U.S. Should Do More to" "Reduce Climate Change (=1)" "' 1.36 `" "International Community Should Do" "More to Reduce Climate Change (=1)" "') xtitle("Hurricane Exposure x Post", size(vlarge)) xscale(titlegap(*5)) xlabel(-.04(.04).2) xmtick(-.04(.02).2) legend(off) yline(.77, lcolor(black) lwidth(*1.5)) xline(0, lcolor(cranberry) lpatt(shortdash))
graph export "${figs}/climchg_action.png", replace
esttab using "${tabs}/ClimChg_Action_SubItems.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{9}{c}{\shortstack{\textbf{Policy Action on Climate Change}}} \\ \cline{3-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Index" "\shortstack{Tax Increase to\\Fund Programs to\\Reduce Climate Change}" "\shortstack{Vote for a Politician\\Who Promised to\\Reduce Climate Change}" "\shortstack{U.S. Should Do\\More to Reduce\\Climate Change}" "\shortstack{Int'l Community Should\\Do More to Reduce\\Climate Change}") extracols(1 2 3 4 5) replace

drop exposure
eststo clear

********************************************************************************
* 					 CLIMATE MITIGATION AND ADAPTATION			   			   *
********************************************************************************
** Table SI-7 of the Additional Supplement

eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe climchg_mitigation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=ian_windswath2
eststo: reghdfe climchg_mitigation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Windswath"

replace exposure=i_ian_surge
eststo: reghdfe climchg_mitigation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Storm Surge"

replace exposure=fema_aid
eststo: reghdfe climchg_mitigation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "FEMA Aid"

replace exposure=hurricane_index1
eststo: reghdfe climchg_adaptation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=ian_windswath2
eststo: reghdfe climchg_adaptation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Windswath"

replace exposure=i_ian_surge
eststo: reghdfe climchg_adaptation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Storm Surge"

replace exposure=fema_aid
eststo: reghdfe climchg_adaptation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "FEMA Aid"

esttab using "${tabs}/Mitigation_Adaptation.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{Climate Change Mitigation Policies}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Climate Change Adaptation Policies}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 5) replace
esttab using "${tabs}/Mitigation_Adaptation_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{4}{c}{\shortstack{\textbf{Climate Change Mitigation Policies}}} & & \multicolumn{4}{c}{\shortstack{\textbf{Climate Change Adaptation Policies}}} \\ \cline{3-6} \cline{8-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 5) replace

drop exposure
eststo clear
** Figure 4e of the Main Paper

gen exposure=hurricane_index1 
eststo: reghdfe climchg_mitigation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_carbon_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_cleanenergy_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_minrestrict_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_fueleff_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est1, msymbol(O) mfcolor(black%0) mlcolor(black%0) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black%0) mlabpos(1) mlabgap(*2) ciopts(lcolor(black%0 black%0) lwidth(.55 1.1))) (est2, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est3, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est4, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est5, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.63 `" "Climate Change Mitigation Policies" "' .93 `" "Carbon Tax (=1)" "' 1.07 `" "Tax to Fund" "Clean Energy (=1)" "' 1.22 `" "Restrictions on Fossil" "Fuel Extraction (=1)" "' 1.36 `" "Stricter Fuel" "Efficiency Standards (=1)" "') xtitle("Hurricane Exposure x Post", size(vlarge)) xscale(titlegap(*5)) xlabel(-.04(.04).2) xmtick(-.04(.02).2) legend(off) yline(.77, lcolor(black) lwidth(*1.5)) xline(0, lcolor(cranberry) lpatt(shortdash))
graph export "${figs}/climchg_mitigation.png", replace
esttab using "${tabs}/Climate_Mitigation_SubItems.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{9}{c}{\shortstack{\textbf{Climate Change Mitigation Policies}}} \\ \cline{3-11} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Index" "\shortstack{Carbon\\Tax}" "\shortstack{Tax to Fund\\Clean Energy}" "\shortstack{Restrictions on Fossil\\Fuel Extraction}" "\shortstack{Stricter Fuel\\Efficiency Standards}") extracols(1 2 3 4 5) replace

eststo clear
** Figure 4f of the Main Paper

eststo: reghdfe climchg_adaptation 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_milclim_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_coast_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_flood_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_streetpump_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe tax_windows_any 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est1, msymbol(O) mfcolor(black%0) mlcolor(black%0) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black%0) mlabpos(1) mlabgap(*2) ciopts(lcolor(black%0 black%0) lwidth(.55 1.1))) (est2, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est3, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est4, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est5, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est6, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.63 `" "Climate Change Adaptation Policies" "' .88 `" "Protecting Military Bases" "From Climate Impacts (=1)" "' 1 `" "Strengthening Coastlines (=1)" "' 1.13 `" "Flood Retrofitting (=1)" "' 1.25 `" "Raising Streets and" "Installing Pumping Stations (=1)" "' 1.38 `" "Requiring Weather-" "Proofed Windows (=1)" "') xtitle("Hurricane Exposure x Post", size(vlarge)) xscale(titlegap(*5)) xlabel(-.04(.04).2) xmtick(-.04(.02).2) legend(off) yline(.755, lcolor(black) lwidth(*1.5)) xline(0, lcolor(cranberry) lpatt(shortdash))
graph export "${figs}/climchg_adaptation.png", replace
esttab using "${tabs}/Climate_Adaptation_SubItems.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{13}{c}} \hline & \\ & & \multicolumn{11}{c}{\shortstack{\textbf{Climate Change Adaptation Policies}}} \\ \cline{3-13} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Index" "\shortstack{Protecting Military Bases\\From Climate Impacts}" "\shortstack{Strengthening\\Coastlines}" "\shortstack{Flood\\Retrofitting}" "\shortstack{Raising Streets and\\Installing Pumping\\Stations}" "\shortstack{Requiring Weather-\\Proofed Windows}") extracols(1 2 3 4 5 6) replace

drop exposure
eststo clear

********************************************************************************
*			 					 CLIMATE SCIENCE				   			   *
********************************************************************************
** Table SI-8 of the Additional Supplement

gen exposure=hurricane_index1
eststo: reghdfe climchg_science 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

replace exposure=ian_windswath2
eststo: reghdfe climchg_science 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Windswath"

replace exposure=i_ian_surge
eststo: reghdfe climchg_science 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Storm Surge"

replace exposure=fema_aid
eststo: reghdfe climchg_science 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "FEMA Aid"

esttab using "${tabs}/Climate_Science.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure) stats(N aic Exposure cov County Date Demographic, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE" "\hspace{3mm}Demographic Covariates")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\shortstack{\textbf{Science of Climate Change}}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 2 3 4) replace
esttab using "${tabs}/Climate_Science_Covariates.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\shortstack{\textbf{Science of Climate Change}}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) nomtit extracols(1 2 3 4) replace

drop exposure
eststo clear
** Figure 4g of the Main Paper

gen exposure=hurricane_index1
eststo: reghdfe climchg_science 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe anthropogenic_agree 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe hurricane_cause_agree 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe hurricane_worse_agree 1.post#c.exposure $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est1, msymbol(O) mfcolor(black%0) mlcolor(black%0) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabsize(large) mlabcolor(black%0) mlabpos(1) mlabgap(*2) ciopts(lcolor(black%0 black%0) lwidth(.55 1.1))) (est2, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est3, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))) (est4, msymbol(O) mfcolor(gs9) mlcolor(gs9) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabsize(large) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(gs9 gs9) lwidth(.55 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.67 `" "Science of Climate Change" "' 1 `" "Human Activities Cause" "Climate Change (=1)" "' 1.17 `" "Climate Change" "Causes Hurricanes (=1)" "' 1.34 `" "Climate Change" "Worsens Hurricanes (=1)" "') xtitle("Hurricane Exposure x Post", size(vlarge)) xscale(titlegap(*5)) xlabel(-.04(.04).2) xmtick(-.04(.02).2) legend(off) yline(.835, lcolor(black) lwidth(*1.5)) xline(0, lcolor(cranberry) lpatt(shortdash))
graph export "${figs}/climate_science.png", replace
esttab using "${tabs}/Climate_Science_SubItems.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\shortstack{\textbf{Science of Climate Change}}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("Index" "\shortstack{Human Activities Cause\\Climate Change}" "\shortstack{Climate Change\\Causes Hurricanes}" "\shortstack{Climate Change\\Worsens Hurricanes}") extracols(1 2 3 4) replace

drop exposure
eststo clear

********************************************************************************
* 								 MIGRATION INTENTIONS						   *
********************************************************************************
** Figure A-4 of the Main Appendix

eststo clear

gen exposure=hurricane_index1 
eststo: reghdfe willmoveclimate 1.post#c.exposure $core $move if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe willmove 1.post#c.exposure $core $move if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe willmovefurther 1.post#c.exposure $core $move if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"
 
coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(12) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est2, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(11) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))) (est3, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(11) mlabgap(*2) ciopts(lcolor(black black) lwidth(.55 1.1))), keep(1.post#c.exposure) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(1 `" "Planning to Move" "in Next 6 Years (=1)" "' 1.25 `" "Planning to Move Further" "from the Coast (=1)" "' .75 `" "Climate Change Increases My" "Likelihood of Moving in the Future (=1)" "') xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xlabel(-.08(.02).1) xmtick(-.08(.01).1) legend(off) xline(0, lcolor(cranberry) lpatt(shortdash))
graph export "${figs}/migration_intention.png", replace
esttab using "${tabs}/Migration_Intention.tex", style(tex) b(3) se(3) nonotes keep(1.post#c.exposure republican democrat woman highschool college age owner timecomm community_links) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.post#c.exposure "& \\ Hurricane Exposure x Post" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age" owner "Home Owner" timecomm "Time Living in Community" community_links "Member of Community Groups") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{7}{c}} \hline & \\ & & \multicolumn{5}{c}{\shortstack{\textbf{Migration Intentions}}} \\ \cline{3-7} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{Climate Change Increases My\\Likelihood of Moving in the Future}" "\shortstack{Planning to Move\\in Next 6 Years}" "\shortstack{Planning to Move Further\\From the Coast}") extracols(1 2 3) replace

drop exposure
eststo clear

********************************************************************************
* 									 EVENT STUDY							   *
********************************************************************************
** Table D-6 of the Additional Supplement

gen sdaysrelative=.
replace sdaysrelative=0 if time2ian==0
replace sdaysrelative=1 if time2ian==1
replace sdaysrelative=2 if time2ian==5
replace sdaysrelative=3 if time2ian==6
replace sdaysrelative=4 if time2ian==9
replace sdaysrelative=5 if time2ian==14
replace sdaysrelative=6 if time2ian==23
replace sdaysrelative=7 if time2ian==29
replace sdaysrelative=-1 if time2ian==-1
replace sdaysrelative=-2 if time2ian==-2
replace sdaysrelative=-3 if time2ian==-9
replace sdaysrelative=-4 if time2ian==-22
replace sdaysrelative=-5 if time2ian==-30
replace sdaysrelative=-5 if time2ian==-31
replace sdaysrelative=-5 if time2ian==-32
replace sdaysrelative=-6 if time2ian==-33
replace sdaysrelative=-6 if time2ian==-35
replace sdaysrelative=-6 if time2ian==-36
replace sdaysrelative=-6 if time2ian==-37
replace sdaysrelative=-7 if time2ian==-44
replace sdaysrelative=-7 if time2ian==-47
replace sdaysrelative=-7 if time2ian==-48
label var sdaysrelative "Survey Rounds Relative to Treatment"
tab sdaysrelative, gen(ll)
gen time=sdaysrelative+8
gen pre_mean=.
gen post_mean=.

eststo clear

eststo: reghdfe climmig_importance 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climmig_action 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_importance 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_action 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_mitigation 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_adaptation 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

eststo: reghdfe climchg_science 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)
estadd local County "Yes"
estadd local Date "Yes"
estadd local Demographic "Yes"
estadd local Exposure "Index"

esttab using "${tabs}/EventStudy_Climate.tex", style(tex) b(3) se(3) nonotes keep(1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 republican democrat woman highschool college age) stats(N aic Exposure cov County Date, fmt(0 3 3) labels("Observations" "AIC" "& \\Exposure Measure:" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}County FE" "\hspace{3mm}Date of Survey FE")) varlabels(1.ll1#c.hurricane_index1 "& \\ Rounds Relative to Treatment: -7" 1.ll2#c.hurricane_index1 "& \\ Rounds Relative to Treatment: -6" 1.ll3#c.hurricane_index1 "& \\ Rounds Relative to Treatment: -5" 1.ll4#c.hurricane_index1 "& \\ Rounds Relative to Treatment: -4" 1.ll5#c.hurricane_index1 "& \\ Rounds Relative to Treatment: -3" 1.ll6#c.hurricane_index1 "& \\ Rounds Relative to Treatment: -2" 1.ll8#c.hurricane_index1 "& \\ Rounds Relative to Treatment: 0" 1.ll9#c.hurricane_index1 "& \\ Rounds Relative to Treatment: 1" 1.ll10#c.hurricane_index1 "& \\ Rounds Relative to Treatment: 2" 1.ll11#c.hurricane_index1 "& \\ Rounds Relative to Treatment: 3" 1.ll12#c.hurricane_index1 "& \\ Rounds Relative to Treatment: 4" 1.ll13#c.hurricane_index1 "& \\ Rounds Relative to Treatment: 5" 1.ll14#c.hurricane_index1 "& \\ Rounds Relative to Treatment: 6" 1.ll15#c.hurricane_index1 "& \\ Rounds Relative to Treatment: 7" republican "& \\ Republican" democrat "Democrat" woman "Woman" highschool "High School Graduate" college "College Graduate" age "Age") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{15}{c}} \hline & \\ & & \multicolumn{13}{c}{\shortstack{\textbf{Focal Indices}}} \\ \cline{3-15} & \\) posthead(\hline) prefoot() postfoot(\hline \end{tabular}) mtit("\shortstack{Issue Importance of\\Climate Migration}" "\shortstack{Policy Action on\\Climate Migration}" "\shortstack{Issue Importance of\\Climate Change}" "\shortstack{Policy Action on\\Climate Change}" "\shortstack{Climate Change\\Mitigation Policies}" "\shortstack{Climate Change\\Adaptation Policies}" "\shortstack{Science of\\Climate Change}") extracols(1 2 3 4 5 6 7) replace

eststo clear

qui reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
gen esample=(e(sample))

gen bjsvar=0
replace bjsvar=. if hurricane_index1<=-.6112757

*** Climate Migration Importance Event Study ***
** Figure SI-1a of the Additional Supplement

eststo clear

eststo: reghdfe climmig_importance 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)

replace pre_mean=.021264
replace post_mean=.1224481

coefplot, vert drop($core _cons) ci(95 90) xline(6.5, lcolor(gs8) lpatt(solid)) yline(0, lcolor(cranberry) lpatt(shortdash)) msymbol(O) msize(vlarge) mcolor(black) ciopts(lwidth(.55 1.1) lcolor(black black)) xlabel(1 "-7" 2 "-6" 3 "-5" 4 "-4" 5 "-3" 6 "-2" 7 "0" 8 "1" 9 "2" 10 "3" 11 "4" 12 "5" 13 "6" 14 "7") xtitle("Survey Rounds Relative to Treatment", size(medlarge)) xscale(titlegap(*5)) ylabel(-.4(.2).6) ymtick(-.4(.1).6) ytitle("Issue Importance of Climate Migration", size(medlarge)) yscale(titlegap(*-1)) nodrop addplot(line pre_mean time if inrange(time,-1,6.5), lcolor(gs7%50) lpatt(solid) || line post_mean time if inrange(time,7,14), lcolor(gs7%50) lpatt(solid))
graph export "${figs}/estudy_climmig_importance.png", replace

eststo clear

eststo: did_imputation climmig_importance geoid sdaysrelative bjsvar if esample==1 [aw=eweight], cluster(geoid) autosample fe(geoid sdaysrelative) minn(0) controls($core) pretrends(6)
event_plot est1, plottype(connected) ciplottype(rarea) graph_opt(title("Climate Migration Issue Importance") xtitle("Time Relative to Treatment") ytitle("Effect of Hurricane Exposure") xlabel(-6(1)-1) ylabel(-1(.5)1) ymtick(-1(.25)1) xline(0, lcolor(gs8) lpattern(dash)) yline(0, lcolor(cranberry) lp(shortdash)) graphregion(color(white)) bgcolor(white) ylabel(, angle(horizontal))) lag_opt(color(navy)) lead_opt(color(black%100) msymbol(O)) lag_ci_opt(color(navy%45 navy%45)) lead_ci_opt(color(gs10%45 gs10%45)) legend_opt(off)
graph export "${figs}/pretrend_climmig_importance.png", replace

eststo clear

*** Climate Migration Policy Event Study ***
** Figure SI-1b of the Additional Supplement

eststo clear

eststo: reghdfe climmig_action 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)

replace pre_mean=.0640315
replace post_mean=.1819173

coefplot, vert drop($core _cons) ci(95 90) xline(6.5, lcolor(gs8) lpatt(solid)) yline(0, lcolor(cranberry) lpatt(shortdash)) msymbol(O) msize(vlarge) mcolor(black) ciopts(lwidth(.55 1.1) lcolor(black black)) xlabel(1 "-7" 2 "-6" 3 "-5" 4 "-4" 5 "-3" 6 "-2" 7 "0" 8 "1" 9 "2" 10 "3" 11 "4" 12 "5" 13 "6" 14 "7") xtitle("Survey Rounds Relative to Treatment", size(medlarge)) xscale(titlegap(*5)) ylabel(-.4(.2).8) ymtick(-.4(.1).8) ytitle("Policy Action on Climate Migration", size(medlarge)) yscale(titlegap(*-1)) nodrop addplot(line pre_mean time if inrange(time,-1,6.5), lcolor(gs7%50) lpatt(solid) || line post_mean time if inrange(time,7,14), lcolor(gs7%50) lpatt(solid))
graph export "${figs}/estudy_climmig_action.png", replace

eststo clear

eststo: did_imputation climmig_action geoid sdaysrelative bjsvar if esample==1 [aw=eweight], cluster(geoid) autosample fe(geoid sdaysrelative) minn(0) controls($core) pretrends(6)
event_plot est1, plottype(connected) ciplottype(rarea) graph_opt(title("Climate Migration Policy Action") xtitle("Time Relative to Treatment") ytitle("Effect of Hurricane Exposure") xlabel(-6(1)-1) ylabel(-1.5(.5)1) ymtick(-1.5(.25)1) xline(0, lcolor(gs8) lpattern(dash)) yline(0, lcolor(cranberry) lp(shortdash)) graphregion(color(white)) bgcolor(white) ylabel(, angle(horizontal))) lag_opt(color(navy)) lead_opt(color(black%100) msymbol(O)) lag_ci_opt(color(navy%45 navy%45)) lead_ci_opt(color(gs10%45 gs10%45)) legend_opt(off)
graph export "${figs}/pretrend_climmig_action.png", replace

eststo clear

*** Climate Change Importance Event Study ***
** Figure SI-1c of the Additional Supplement

eststo clear

eststo: reghdfe climchg_importance 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)

replace pre_mean=.1435088
replace post_mean=.2640303

coefplot, vert drop($core _cons) ci(95 90) xline(6.5, lcolor(gs8) lpatt(solid)) yline(0, lcolor(cranberry) lpatt(shortdash)) msymbol(O) msize(vlarge) mcolor(black) ciopts(lwidth(.55 1.1) lcolor(black black)) xlabel(1 "-7" 2 "-6" 3 "-5" 4 "-4" 5 "-3" 6 "-2" 7 "0" 8 "1" 9 "2" 10 "3" 11 "4" 12 "5" 13 "6" 14 "7") xtitle("Survey Rounds Relative to Treatment", size(medlarge)) xscale(titlegap(*5)) ylabel(-.4(.2).8) ymtick(-.4(.1).8) ytitle("Issue Importance of Climate Change", size(medlarge)) yscale(titlegap(*-1)) nodrop addplot(line pre_mean time if inrange(time,-1,6.5), lcolor(gs7%50) lpatt(solid) || line post_mean time if inrange(time,7,14), lcolor(gs7%50) lpatt(solid))
graph export "${figs}/estudy_climchg_importance.png", replace

eststo clear

eststo: did_imputation climchg_importance geoid sdaysrelative bjsvar if esample==1 [aw=eweight], cluster(geoid) autosample fe(geoid sdaysrelative) minn(0) controls($core) pretrends(6)
event_plot est1, plottype(connected) ciplottype(rarea) graph_opt(title("Climate Change Issue Importance") xtitle("Time Relative to Treatment") ytitle("Effect of Hurricane Exposure") xlabel(-6(1)-1) ylabel(-1.5(.5)1) ymtick(-1.5(.25)1) xline(0, lcolor(gs8) lpattern(dash)) yline(0, lcolor(cranberry) lp(shortdash)) graphregion(color(white)) bgcolor(white) ylabel(, angle(horizontal))) lag_opt(color(navy)) lead_opt(color(black%100) msymbol(O)) lag_ci_opt(color(navy%45 navy%45)) lead_ci_opt(color(gs10%45 gs10%45)) legend_opt(off)
graph export "${figs}/pretrend_climchg_importance.png", replace

eststo clear

*** Climate Change Policy Event Study ***
** Figure SI-1d of the Additional Supplement

eststo clear

eststo: reghdfe climchg_action 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)

replace pre_mean=.0104087
replace post_mean=.1321758

coefplot, vert drop($core _cons) ci(95 90) xline(6.5, lcolor(gs8) lpatt(solid)) yline(0, lcolor(cranberry) lpatt(shortdash)) msymbol(O) msize(vlarge) mcolor(black) ciopts(lwidth(.55 1.1) lcolor(black black)) xlabel(1 "-7" 2 "-6" 3 "-5" 4 "-4" 5 "-3" 6 "-2" 7 "0" 8 "1" 9 "2" 10 "3" 11 "4" 12 "5" 13 "6" 14 "7") xtitle("Survey Rounds Relative to Treatment", size(medlarge)) xscale(titlegap(*5)) ylabel(-.6(.2).8) ymtick(-.6(.1).8) ytitle("Policy Action on Climate Change", size(medlarge)) yscale(titlegap(*-30)) nodrop addplot(line pre_mean time if inrange(time,-1,6.5), lcolor(gs7%50) lpatt(solid) || line post_mean time if inrange(time,7,14), lcolor(gs7%50) lpatt(solid))
graph export "${figs}/estudy_climchg_action.png", replace

eststo clear

eststo: did_imputation climchg_action geoid sdaysrelative bjsvar if esample==1 [aw=eweight], cluster(geoid) autosample fe(geoid sdaysrelative) minn(0) controls($core) pretrends(6)
event_plot est1, plottype(connected) ciplottype(rarea) graph_opt(title("Climate Change Policy Action") xtitle("Time Relative to Treatment") ytitle("Effect of Hurricane Exposure") xlabel(-6(1)-1) ylabel(-1.5(.5)1) ymtick(-1.5(.25)1) xline(0, lcolor(gs8) lpattern(dash)) yline(0, lcolor(cranberry) lp(shortdash)) graphregion(color(white)) bgcolor(white) ylabel(, angle(horizontal))) lag_opt(color(navy)) lead_opt(color(black%100) msymbol(O)) lag_ci_opt(color(navy%45 navy%45)) lead_ci_opt(color(gs10%45 gs10%45)) legend_opt(off)
graph export "${figs}/pretrend_climchg_action.png", replace

eststo clear

*** Climate Change Mitigation Event Study ***
** Figure SI-2a of the Additional Supplement

eststo clear

eststo: reghdfe climchg_mitigation 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)

replace pre_mean=.076355
replace post_mean=.1889062

coefplot, vert drop($core _cons) ci(95 90) xline(6.5, lcolor(gs8) lpatt(solid)) yline(0, lcolor(cranberry) lpatt(shortdash)) msymbol(O) msize(vlarge) mcolor(black) ciopts(lwidth(.55 1.1) lcolor(black black)) xlabel(1 "-7" 2 "-6" 3 "-5" 4 "-4" 5 "-3" 6 "-2" 7 "0" 8 "1" 9 "2" 10 "3" 11 "4" 12 "5" 13 "6" 14 "7") xtitle("Survey Rounds Relative to Treatment", size(medlarge)) xscale(titlegap(*5)) ylabel(-.6(.2).8) ymtick(-.6(.1).8) ytitle("Climate Change Mitigation Policies", size(medlarge)) yscale(titlegap(*-30)) nodrop addplot(line pre_mean time if inrange(time,-1,6.5), lcolor(gs7%50) lpatt(solid) || line post_mean time if inrange(time,7,14), lcolor(gs7%50) lpatt(solid))
graph export "${figs}/estudy_climchg_mitigation.png", replace

eststo clear

eststo: did_imputation climchg_mitigation geoid sdaysrelative bjsvar if esample==1 [aw=eweight], cluster(geoid) autosample fe(geoid sdaysrelative) minn(0) controls($core) pretrends(6)
event_plot est1, plottype(connected) ciplottype(rarea) graph_opt(title("Climate Change Mitigation Policies") xtitle("Time Relative to Treatment") ytitle("Effect of Hurricane Exposure") xlabel(-6(1)-1) ylabel(-1(.5)1.5) ymtick(-1(.25)1.5) xline(0, lcolor(gs8) lpattern(dash)) yline(0, lcolor(cranberry) lp(shortdash)) graphregion(color(white)) bgcolor(white) ylabel(, angle(horizontal))) lag_opt(color(navy)) lead_opt(color(black%100) msymbol(O)) lag_ci_opt(color(navy%45 navy%45)) lead_ci_opt(color(gs10%45 gs10%45)) legend_opt(off)
graph export "${figs}/pretrend_climchg_mitigation.png", replace

eststo clear

*** Climate Change Adaptation Event Study ***
** Figure SI-2b of the Additional Supplement

eststo clear

eststo: reghdfe climchg_adaptation 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)

replace pre_mean=-.000798
replace post_mean=.1469445

coefplot, vert drop($core _cons) ci(95 90) xline(6.5, lcolor(gs8) lpatt(solid)) yline(0, lcolor(cranberry) lpatt(shortdash)) msymbol(O) msize(vlarge) mcolor(black) ciopts(lwidth(.55 1.1) lcolor(black black)) xlabel(1 "-7" 2 "-6" 3 "-5" 4 "-4" 5 "-3" 6 "-2" 7 "0" 8 "1" 9 "2" 10 "3" 11 "4" 12 "5" 13 "6" 14 "7") xtitle("Survey Rounds Relative to Treatment", size(medlarge)) xscale(titlegap(*5)) ylabel(-.6(.2).8) ymtick(-.6(.1).8) ytitle("Climate Change Adaptation Policies", size(medlarge)) yscale(titlegap(*-30)) nodrop addplot(line pre_mean time if inrange(time,-1,6.5), lcolor(gs7%50) lpatt(solid) || line post_mean time if inrange(time,7,14), lcolor(gs7%50) lpatt(solid))
graph export "${figs}/estudy_climchg_adaptation.png", replace

eststo clear

eststo: did_imputation climchg_adaptation geoid sdaysrelative bjsvar if esample==1 [aw=eweight], cluster(geoid) autosample fe(geoid sdaysrelative) minn(0) controls($core) pretrends(6)
event_plot est1, plottype(connected) ciplottype(rarea) graph_opt(title("Climate Change Adaptation Policies") xtitle("Time Relative to Treatment") ytitle("Effect of Hurricane Exposure") xlabel(-6(1)-1) ylabel(-1(.5)1.5) ymtick(-1(.25)1.5) xline(0, lcolor(gs8) lpattern(dash)) yline(0, lcolor(cranberry) lp(shortdash)) graphregion(color(white)) bgcolor(white) ylabel(, angle(horizontal))) lag_opt(color(navy)) lead_opt(color(black%100) msymbol(O)) lag_ci_opt(color(navy%45 navy%45)) lead_ci_opt(color(gs10%45 gs10%45)) legend_opt(off)
graph export "${figs}/pretrend_climchg_adaptation.png", replace

eststo clear

*** Climate Change Science Event Study ***
** Figure SI-2C of the Additional Supplement

eststo clear

eststo: reghdfe climchg_science 1.ll1#c.hurricane_index1 1.ll2#c.hurricane_index1 1.ll3#c.hurricane_index1 1.ll4#c.hurricane_index1 1.ll5#c.hurricane_index1 1.ll6#c.hurricane_index1 1.ll8#c.hurricane_index1 1.ll9#c.hurricane_index1 1.ll10#c.hurricane_index1 1.ll11#c.hurricane_index1 1.ll12#c.hurricane_index1 1.ll13#c.hurricane_index1 1.ll14#c.hurricane_index1 1.ll15#c.hurricane_index1 $core if estsample==1 & inrange(sdaysrelative,-7,15) [pw=eweight], cluster(geoid) abs(geoid start_date)

replace pre_mean=.0665257
replace post_mean=.1872109

coefplot, vert drop($core _cons) ci(95 90) xline(6.5, lcolor(gs8) lpatt(solid)) yline(0, lcolor(cranberry) lpatt(shortdash)) msymbol(O) msize(vlarge) mcolor(black) ciopts(lwidth(.55 1.1) lcolor(black black)) xlabel(1 "-7" 2 "-6" 3 "-5" 4 "-4" 5 "-3" 6 "-2" 7 "0" 8 "1" 9 "2" 10 "3" 11 "4" 12 "5" 13 "6" 14 "7") xtitle("Survey Rounds Relative to Treatment", size(medlarge)) xscale(titlegap(*5)) ylabel(-.6(.2).8) ymtick(-.6(.1).8) ytitle("Science of Climate Change", size(medlarge)) yscale(titlegap(*-30)) nodrop addplot(line pre_mean time if inrange(time,-1,6.5), lcolor(gs7%50) lpatt(solid) || line post_mean time if inrange(time,7,14), lcolor(gs7%50) lpatt(solid))
graph export "${figs}/estudy_climchg_science.png", replace

eststo clear

eststo: did_imputation climchg_science geoid sdaysrelative bjsvar if esample==1 [aw=eweight], cluster(geoid) autosample fe(geoid sdaysrelative) minn(0) controls($core) pretrends(6)
event_plot est1, plottype(connected) ciplottype(rarea) graph_opt(title("Science of Climate Change") xtitle("Time Relative to Treatment") ytitle("Effect of Hurricane Exposure") xlabel(-6(1)-1) ylabel(-1.5(.5)1) ymtick(-1.5(.25)1) xline(0, lcolor(gs8) lpattern(dash)) yline(0, lcolor(cranberry) lp(shortdash)) graphregion(color(white)) bgcolor(white) ylabel(, angle(horizontal))) lag_opt(color(navy)) lead_opt(color(black%100) msymbol(O)) lag_ci_opt(color(navy%45 navy%45)) lead_ci_opt(color(gs10%45 gs10%45)) legend_opt(off)
graph export "${figs}/pretrend_climchg_science.png", replace

eststo clear
drop esample

qui reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
gen esample=(e(sample))

foreach x of var climmig_importance climmig_action climchg_importance climchg_action climchg_mitigation climchg_adaptation climchg_science {

matrix storeDID`x' = J(7, 4, .) 
matrix list storeDID`x'  
forval i=1/7{
gen dateDID=.
replace dateDID=0 if time==`i'+1
replace dateDID=1 if time==`i'
diff `x' if esample==1 [aw=eweight], t(binary_index1) p(dateDID) cov($core)
matrix storeDID`x'[`i', 1] = `i'
matrix storeDID`x'[`i', 2] = `r(did)'
matrix storeDID`x'[`i', 3] = `r(se_dd)'
matrix storeDID`x'[`i', 4] = `r(did)'/`r(se_dd)'
drop dateDID
}
matrix list storeDID`x'  
}
drop _diff

********************************************************************************
* 							 DROPPING POST PERIODS							   *
********************************************************************************
** Result not presented in paper

eststo clear

eststo: reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=0 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=2 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=3 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=4 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=5 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=6 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=7 [pw=eweight], cluster(geoid) abs(geoid start_date)

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est2, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est3, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est4, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est5, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est6, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est7, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est8, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est9, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))), keep(1.post#c.hurricane_index1) legend(off) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.6 "Benchmark" .7 "Drop {it:t + 0}" .8 "Drop {it:t + 1}" .9 "Drop {it:t + 2}" 1 "Drop {it:t + 3}" 1.1 "Drop {it:t + 4}" 1.2 "Drop {it:t + 5}" 1.3 "Drop {it:t + 6}" 1.4 "Drop {it:t + 7}") title("Climate Migration Issue Importance") xlabel(0(.05).25) xmtick(0(.01).25)
graph export "${figs}/drop_climmig_importance.png", replace

eststo clear

eststo: reghdfe climmig_action 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=0 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=2 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=3 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=4 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=5 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=6 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climmig_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=7 [pw=eweight], cluster(geoid) abs(geoid start_date)

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est2, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est3, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est4, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est5, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est6, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est7, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est8, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est9, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))), keep(1.post#c.hurricane_index1) legend(off) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.6 "Benchmark" .7 "Drop {it:t + 0}" .8 "Drop {it:t + 1}" .9 "Drop {it:t + 2}" 1 "Drop {it:t + 3}" 1.1 "Drop {it:t + 4}" 1.2 "Drop {it:t + 5}" 1.3 "Drop {it:t + 6}" 1.4 "Drop {it:t + 7}") title("Climate Migration Policy Action") xlabel(0(.05).25) xmtick(0(.01).25)
graph export "${figs}/drop_climmig_action.png", replace

eststo clear

eststo: reghdfe climchg_importance 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=0 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=2 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=3 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=4 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=5 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=6 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_importance 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=7 [pw=eweight], cluster(geoid) abs(geoid start_date)

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est2, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est3, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est4, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est5, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est6, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est7, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est8, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est9, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))), keep(1.post#c.hurricane_index1) legend(off) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.6 "Benchmark" .7 "Drop {it:t + 0}" .8 "Drop {it:t + 1}" .9 "Drop {it:t + 2}" 1 "Drop {it:t + 3}" 1.1 "Drop {it:t + 4}" 1.2 "Drop {it:t + 5}" 1.3 "Drop {it:t + 6}" 1.4 "Drop {it:t + 7}") title("Climate Change Issue Importance") xlabel(0(.05).25) xmtick(0(.01).25)
graph export "${figs}/drop_climchg_importance.png", replace

eststo clear

eststo: reghdfe climchg_action 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=0 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=2 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=3 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=4 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=5 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=6 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_action 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=7 [pw=eweight], cluster(geoid) abs(geoid start_date)

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est2, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est3, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est4, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est5, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est6, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est7, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est8, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est9, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))), keep(1.post#c.hurricane_index1) legend(off) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.6 "Benchmark" .7 "Drop {it:t + 0}" .8 "Drop {it:t + 1}" .9 "Drop {it:t + 2}" 1 "Drop {it:t + 3}" 1.1 "Drop {it:t + 4}" 1.2 "Drop {it:t + 5}" 1.3 "Drop {it:t + 6}" 1.4 "Drop {it:t + 7}") title("Climate Change Policy Action") xlabel(0(.05).25) xmtick(0(.01).25)
graph export "${figs}/drop_climchg_action.png", replace

eststo clear

eststo: reghdfe climchg_mitigation 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_mitigation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=0 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_mitigation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_mitigation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=2 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_mitigation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=3 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_mitigation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=4 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_mitigation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=5 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_mitigation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=6 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_mitigation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=7 [pw=eweight], cluster(geoid) abs(geoid start_date)

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est2, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est3, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est4, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est5, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est6, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est7, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est8, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est9, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))), keep(1.post#c.hurricane_index1) legend(off) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.6 "Benchmark" .7 "Drop {it:t + 0}" .8 "Drop {it:t + 1}" .9 "Drop {it:t + 2}" 1 "Drop {it:t + 3}" 1.1 "Drop {it:t + 4}" 1.2 "Drop {it:t + 5}" 1.3 "Drop {it:t + 6}" 1.4 "Drop {it:t + 7}") title("Climate Change Mitigation Policies") xlabel(-.05(.05).25) xmtick(-0.05(.01).25)
graph export "${figs}/drop_climchg_mitigation.png", replace

eststo clear

eststo: reghdfe climchg_adaptation 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_adaptation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=0 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_adaptation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_adaptation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=2 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_adaptation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=3 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_adaptation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=4 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_adaptation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=5 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_adaptation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=6 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_adaptation 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=7 [pw=eweight], cluster(geoid) abs(geoid start_date)

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est2, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est3, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est4, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est5, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est6, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est7, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est8, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est9, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel("{it:p} = " + string(@pval,"%9.3f")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))), keep(1.post#c.hurricane_index1) legend(off) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.6 "Benchmark" .7 "Drop {it:t + 0}" .8 "Drop {it:t + 1}" .9 "Drop {it:t + 2}" 1 "Drop {it:t + 3}" 1.1 "Drop {it:t + 4}" 1.2 "Drop {it:t + 5}" 1.3 "Drop {it:t + 6}" 1.4 "Drop {it:t + 7}") title("Climate Change Adaptation Policies") xlabel(-.05(.05).25) xmtick(-0.05(.01).25)
graph export "${figs}/drop_climchg_adaptation.png", replace

eststo clear

eststo: reghdfe climchg_science 1.post#c.hurricane_index1 $core if estsample==1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_science 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=0 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_science 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=1 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_science 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=2 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_science 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=3 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_science 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=4 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_science 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=5 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_science 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=6 [pw=eweight], cluster(geoid) abs(geoid start_date)
eststo: reghdfe climchg_science 1.post#c.hurricane_index1 $core if estsample==1 & sdaysrelative!=7 [pw=eweight], cluster(geoid) abs(geoid start_date)

coefplot (est1, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est2, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est3, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est4, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est5, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est6, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est7, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est8, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))) (est9, msymbol(O) mfcolor(black) mlcolor(black) msize(vlarge) mlabel(cond(@pval<.001, "{it:p} < 0.001", "")) mlabcolor(black) mlabpos(1) mlabgap(*2) ciopts(lcolor(black black) lwidth(.45 1.1))), keep(1.post#c.hurricane_index1) legend(off) ci(95 90) yline(0, lcolor(cranberry) lpatt(shortdash)) xtitle("Hurricane Exposure x Post", size(medium)) xscale(titlegap(*5)) xline(0, lcolor(cranberry) lpatt(shortdash)) ylabel(.6 "Benchmark" .7 "Drop {it:t + 0}" .8 "Drop {it:t + 1}" .9 "Drop {it:t + 2}" 1 "Drop {it:t + 3}" 1.1 "Drop {it:t + 4}" 1.2 "Drop {it:t + 5}" 1.3 "Drop {it:t + 6}" 1.4 "Drop {it:t + 7}") title("Science of Climate Change") xlabel(0(.05).25) xmtick(0(.01).25)
graph export "${figs}/drop_climchg_science.png", replace

drop sdaysrelative time pre_mean post_mean 
eststo clear

********************************************************************************
